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Abstract The German engineer H.A. Janssen gave one of 
the first accounts of the often peculiar behavior of granular 
material in a paper published in German in 1895. From sim- 
ple experiments with corn he inferred the saturation of pres- 
sure with height in a granular system. Subsequently, Janssen 
derived the equivalent of the barometric formula for granu- 
lar material from the main assumption that the walls carry 
part of the weight. The following is a translation of this 
article. The wording is chosen as close as possible to the 
original. While drawings are copied from the original, fig- 
ures displaying data are redone for better readability. The 
translation is complemented by some bibliographical notes 
and an assessment of earlier work, wherein Hagen predicted 
the saturation of pressure with depth in 1852, and Huber- 
Burnand demonstrated that saturation qualitatively as early 
■ as in 1829. We conclude with a brief discussion of more re- 
cent developments resting on Janssen's work. 



1 HA. Janssen, Engineer in Bremen, Germany, 
Experiments on Corn Pressure in Silo Cells, 31st August 
1895 [1] 

During the last decades the export of corn from the corn 
countries of the world to the culture states of Europe has 
made extraordinary progress. Considering the length of the 
route of transport, this development could only have been 
achieved by low production costs, be it due to low wages, be 
it due to extensive application of mechanical equipment both 
for tillage and for harvesting. Further, consider improvements 
of corn depositories and the means of corn transportation. 
Predominantly in the United States of North America stor- 
age silos have proven themselves to be beneficial. Constructed 
on a large scale, these silos are designed to receive the corn 
arriving in railway carriages, to store it, and to dispense it to 
river and sea vessels [2]. The extensive mechanical equip- 
ment for stocking the fruit, weighing it, and loading it onto 
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ships is highly efficient. Some realizations are capable of 
handling several hundreds of tons per hour. The interior of 
such a storage silo - called an elevator in North America 
- consists of a number of vertical tubular cavities serving 
to receive the corn from street-level to the attic. The profile 
of these tubes or cells has been exclusively rectangular with 
common walls from wooden shafts and bottom casing. Re- 
cently the walls are also produced from iron-strengthened 
brickwork, and in six-fold profiles like in honeycombs [2]. 
The grain is introduced through the upper end of the cell 
through a hatch. For discharge, stock transfer, or embarka- 
tion, the bottoms of the cells provide close-able openings. 
The dimensions of the cells vary considerably; the largest 
ones have a capacity of up to around 250 t, in which case 
the height of the corn reaches around 25 m. It is obvious 
that the enormous content of the cell must exert consider- 
able pressure on the side walls and the bottom. While vari- 
ous construction guidebooks indicate the necessity for suf- 
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ficient anchoring of the cell walls, nowhere is anything in 
particular given for the determination of the corn pressure 
against the side walls. The formulas valid for the pressure 
in liquids are not applicable for corn as the friction among 
the grains influences the pressure transmission considerably. 
Likewise, the formulas customary for the calculation of the 
pressure against feed walls are alone for the reason not appli- 
cable for the present case, because the latter are derived for 
straight walls, while the content of a silo cell is surrounded 
by vertical walls from all sides. To my knowledge, the only 
publication about experiments to determine the corn pres- 
sure in silo cells were done in the Engineering on October 
27th, 1882 [3]. These experiments establish that the pressure 
of corn in a silo cell against the support grows with increas- 
ing depth up to a certain value, but does not grow further 
at larger depth. The largest emerging pressure on the area 
depends on the profile of the cell, and is supposedly propor- 
tional to the diameter of the circles inscribed in the profile. 

Since it appeared valuable to me to gain insight about 
the amount of the pressure increase in the silo cell from the 
surface down to arbitrary depth, I performed detailed exper- 
iments, whose results shall be reproduced here. 



Table 1 Runs 1-3 



run 1, cell 1 


run 2, cell 2 


run 3, 


cell 2 


area 20 


•20 cm 
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• 30 cm 
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• 30 cm 
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50.0 
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Four wooden sample cells were produced with quadratic 
profile and side lengths of 20, 30, 40, and 60 cm. During 
the time of the experiment the respective cell is mounted on 
four screws S, Fig. 1 and 2. The lower end of the tube is 
secluded by a well-fitting movable base. This bottom of the 
cell rests on a decimal balance that is counterbalanced by 
weights on the plate G before the experiment. In this cen- 
ter position a foothold is placed underneath the plate G and 
G is loaded with an additional piece of weight. The corn is 
weighed in hollow vessels before being poured into the ap- 
paratus. Now the cell is filled up to the point where plate 
G with the applied weights starts moving upwards from its 
center position. Weighing the remaining corn in the filling 
container, one can determine accurately the content in the 
cell that is causing the now known pressure at the bottom. 
After placing an additional weight on the plate G, the ap- 
paratus is lifted with the screws S to the point where the 
balance is in equilibrium again. Now the experiment is con- 
tinued. In this way it was possible to determine the bottom 
pressure for various filling heights while filling the cell only 
once. 

The results are displayed graphically in Figs. 4 to 7. The 
smooth path of the observed pressure curve for the runs 2 to 
6 lets us conclude that substantial observational errors have 
not occurred. Run number 1 - which was conducted first 
- displays relatively large observational errors, because the 
observers became more proficient in the handling of the in- 
struments only during the course of the experiments. 

Were the content of the cell to consist of a liquid, the 
bottom pressure would be equal to the weight of the cell's 
content. The experiments with corn, however, yield a much 
smaller bottom pressure. This fact can be related to the fric- 
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Table 2 Runs 4-6 



run 4, cell 3 


run 5, 


cell 4 


run 6 


cell 4 


area 40 


■40 cm 


area 60 


■ 60 cm 


area 60 


•60 cm 


content: wheat 


content: wheat 


content: wheat 


spec, weight y = 0.8 


spec, weight y = 0.8 


spec, weight y = 0.8 


bottom 


corn 
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corn 


bottom 
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pressure 


content 


pressure 


content 


pressure 


content 


kg 


kg 


kg 


kg 


kg 


kg 
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52.5 
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66.0 
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85.0 
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Fig. 4 Run 1. 
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Fig. 7 Runs 5 and 6. 
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Fig. 5 Runs 2 and 3. 



so 



S 
o 

o 



100 



200 

corn amount [kg 



Fig. 6 Run 4. 
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tion between corn and cell walls. This friction becomes so 
large at increasing depth that a pressure increase is no longer 
noticeable. Hence, the friction between corn and cell walls 
needs to be equal to the weight of the enclosed layers of 
corn. The magnitude of the corn pressure exerted on the 
encompassing walls in this case - the pressure apparently 
reaches its highest value - can be calculated in the follow- 
ing way. Let p Simdx = largest pressure of the corn against the 
wall, / = friction coefficient between corn and cell walls, 
^ = side of the quadratic cell profile, dh = height of a corn 
layer, y = specific weight of the corn, then Fig. 8 yields the 
equation 



Ps,rmaf4sdh=JS dh , 



JS 

Ps,m-ix = ■ 



(1) 



Assuming a constant coefficient of friction / for variable 
pressures, the maximum pressure exerted on the cell walls by 
the corn for various cell widths of quadratic profile is pro- 
portional to the side length of the cell profile. To determine 
the value of /, the following experiments were performed: 

A flat container, Fig. 9, was filled with wheat and cov- 
ered with a wooden plate. The plate could be pulled across 
the corn with a string that was connected to a spring balance 
F. By putting weights on top the pressure between plate and 
corn could be regulated precisely. The result of these exper- 
iments is given in the following tables. 

The investigations No 10 and 11 were done in the same 
way on different days. The difference of the friction coef- 
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Table 3 Runs 7-9 



run 7 


run 8 


run 9 


pressure 


measured 


pressure 


measured 


pressure 


measured 


lid-corn 


drag 


lid-corn 


drag 


lid-corn 


drag 


kg 


kg 


kg 


kg 


kg 


kg 


5.9 


1.93 


10.9 


3.73 


20.9 


6.98 


5.9 


1.93 


10.9 


3.73 


20.9 


6.73 


5.9 


1.93 


10.9 


3.78 


20.9 


6.78 


5.9 


1.93 


10.9 


3.78 


20.9 


7.23 


5.9 


1.98 


10.9 


3.83 


20.9 


7.28 Fi 


5.9 


1.88 


10.9 


3.73 


20.9 


7.23 


5.9 


2.03 


10.9 


3.78 


20.9 


7.18 


5.9 


2.08 


10.9 


3.78 


20.9 


6.98 


5.9 


1.93 


10.9 


3.73 


20.9 


6.98 i t 


5.9 


1.98 


10.9 


3.73 


20.9 


6.73 


mean drag = 1.96 


mean drag = 3.77 


mean drag = 7.01 


f=1.96/5.9=0.332 


f=3.77/10.9=0.346 


f=7.01/20.9=0.335 




F{p + dp — p) = yFdx — fp s udx 
dp = ydx 



p-dx 



dp 



Table 4 Runs 10-12 



= dx 



run 10 


run 11 


run 12 


pressure 


measured 




pressure 


measured 


pressure 


measured 


lid-corn 


drag 




lid-corn 


drag 


lid-corn 


drag 


kg 


kg 




kg 


kg 


kg 


kg 


30.9 


9.23 




30.9 


10.22 


40.9 


12.53 


30.9 


9.03 




30.9 


10.18 


40.9 


12.48 


30.9 


9.58 




30.9 


10.08 


40.9 


12.38 


30.9 


9.23 




30.9 


9.93 


40.9 


11.93 


30.9 


9.63 




30.9 


9.93 


40.9 


12.28 


30.9 


9.28 




30.9 


10.03 


40.9 


12.48 


30.9 


9.43 




30.9 


10.38 


40.9 


12.28 


30.9 


9.18 




30.9 


10.03 


40.9 


12.53 


mean drag = 9.32 




mean drag = 10.10 


mean drag = 12.36 


f=9.32/30.9=0.302 




f=10. 10/30.9=0.327 


f=12.36/40.9=0.302 
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ficients determined this way - there is a deviation of 8% 
- should be ascribed to different humidity. The maximum 
value for / resulted from run 8 and was 0.346, the minimum 
value from runs 10 and 12 was 0.302. 

These investigations where the pressure was gradually 
increased by a factor of seven show that the value for the 
friction coefficient is not significantly influenced thereby; 
therefore, / can be assumed constant with sufficient accu- 
racy for different filling heights. Provided that the horizontal 
pressure of the corn is proportional to the vertical pressure, 
the latter can be calculated in the following way: 

Let P = total pressure of the corn on the cell floor, p = 
vertical pressure of the corn, p s = horizontal pressure of the 
corn, / = friction coefficient between corn and cell wall, 
K = — ; pK = p s f, s = side length of the quadratic cell 
profile, u = circumference of the cell =4s, F = area of the 
profile = s 2 , x = filling height of the corn in the cell, y = spe- 
cific weight of the content, e = base of the natural logarithm; 



1 AK V 

P=^(l- 
4K V 



-AK 



(2) 



(3) 



In the latter two equations only K is unknown, which 
can be determined from the experiments 1 to 6. In run 4 the 
largest total bottom pressure was found at 63 kg or 3.94 kp 
per 1 qdm. The friction at the circumference in this case for 
a layer of 1 dm is 16y = 16 • 0.8 = 12.8 kg = p ju; 



Psf = 



12.8 12.8 



16 



= 0.8; K = 



P.J 0.8 



3.94 



= 0.203 . 



For simplification I set K = 0.2; y = 0.8. Then Eqs. (2) 
and (3) assume the following form 



P 
P 



(1 



.,-0.8? 



) 



1 



„-0.8* 



(2a) 
(3a) 



Figure 1 1 displays the results of run 4 and the calculated 
corn pressures according to Eq. (3a). Better agreement could 
have been achieved by assuming a larger value of K, consid- 
ering the pressure measurements for smaller filling heights, 
since it can not be ruled out that at the observation of the 
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Fig. 11 Data from run 4 and fit by Eq. (3a) 

largest pressure an observational error has occurred. In the 
same way as described above, K is determined for the other 
runs, hence: 



run 1 


K = 


0.211 


run 2, 3 


K = 


0.235 


run 4 


K = 


0.203 


run 6 


K = 


0.227 



The deviations can be explained by a slight variability 
of /. In runs 2 and 3, / should have been larger by 15.5% 
compared to run 4. The experiments for the determination 
of the value of the friction coefficient result in (according to 
runs 7 to 12) / = 0.302 to 0.346, or 14±% deviation. For 
the calculation of the side pressure against the cell wall one 
can therefore assume that the value of / was around 0.346 
in run 2 and 3 and around 0.302 in run 4. One gets from 
p s = for run 2 and 3: P s = £§f§/? = 0.68/?, and for run 4: 

p s = jp§^p = 0.675/?, or rounded p s = O.lp 

Hereby the mean pressure against a side wall of a cell 
of quadratic profile is determined. However, we can safely 
assume that near the corners of the profile the pressure is 
lower than this mean value, while it is higher in the center 
part of the walls. It is shown in Fig. 12 how the pressure is 
transmitted to the side walls. Thereby it was assumed that 
the distribution of the corn pressure originates radially from 
the center of the cell and exerts a pressure against the wall 
of psina, with a being the angle under which the pressure 
ray hits the wall. The maximum pressure against the wall is 
reached in the middle of the cell wall with around 1.15 of the 
mean pressure, or = 1.15- 0.7 • p w 0.8/?. For the calculation 
of the wall thickness the ansatz of uniform loading can be 
made with sufficient accuracy 

Ps = 0.75/? = ^ • 0.75 (1 - e-°-^) . (4a) 

Even if the above experimental and calculated results are 
sufficient to examine silo bottoms and silo walls for their 
necessary stiffness, a simpler way of calculation might be 
desirable. We are on the safe side for the usual kinds of corn 
if we assume wheat as cell content with a specific weight of 
0.80, set K = 0.20, and calculate the corn pressure at the bot- 
tom according to formulas (2a) and (3a). From the graph in 
Fig. 13 one can obtain the value for the brackets in Eqs. (2a) 
and (3a) = (l - e" 08 ;) = Z for arbitrary f . 




Fig. 13 Calculation of the bottom pressure. 




Fig. 14 Calculation of the side pressures 

From the just determined vertical pressures one could 
get the respective side pressures - which are decisive for the 
dimensioning of the cell walls - by multiplying by 0.75; for 
greater convenience the value of w = 0.75 (l — <?~ 08 «) can 
be obtained from Fig. 14 for arbitrary 

Originally it was the intention of the author to determine 
the side pressure of the corn directly in the experiments. 
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Table 5 Example 1 . Calculation of the bottom pressure for a cell of 
4 x 4 m base: s = 4. Z from Fig. 13, p = sZ, P = s 3 Z, 



filling x/m 


x/s 


Z 


P 


P 


1 


0.25 


0.18 


0.72 


11.5 


2 


0.50 


0.33 


1.32 


21.1 


3 


0.75 


0.45 


1.80 


28.8 


4 


1.00 


0.55 


2.20 


35.2 


5 


1.25 


0.63 


2.52 


40.3 


6 


1.50 


0.70 


2.80 


44.7 


8 


2.0 


0.80 


3.20 


51.2 


10 


2.5 


0.86 


3.44 


55.0 


12 


3.0 


0.91 


3.64 


58.2 


14 


3.5 


0.94 


3.76 


60.1 



Table 6 Example 2. Calculation of the side pressure of a silo cell of 
3 x 3 m base: s = 3. w after Fig. 14, p s = sw, 



filling x[m] 


x/s 


w 


p s [t/m 2 ] 


1.5 


0.50 


0.245 


0.735 


3.0 


1.00 


0.410 


1.23 


4.5 


1.50 


0.530 


1.59 


6.0 


2.00 


0.595 


1.79 


8.0 


2.67 


0.655 


1.97 


10.0 


3.33 


0.690 


2.07 


12.0 


4.00 


0.720 


2.16 



To this end, one of the sample cells was equipped with a 
side lid that was pressed on the cell by suitable weights 
by an angle lever. The pressure necessary to open the lid 
could then be calculated. But while conducting the experi- 
ment the opening of the lid happened so slowly that accu- 
rate results could not be obtained. The same occurred with 
the attempt to measure the bottom pressure at various points 
on the bottom via small lids. In my opinion this feature can 
be traced back to the fact that arching occurs as the profile 
of the corn column narrows, which significantly influences 
the pressure transmission for smaller profiles. The attempt 
to equip the entire experiment from top to bottom with a 
movable cell wall was discarded because substantial incon- 
veniences would have been incurred. Also the calculation 
method described above led to the same goal. 

Besides wheat, further investigations were conducted with 
rye and maize, whose results are in agreement with the ex- 
periments described above. With rye - with a specific weight 
of 0.75 - the pressures fell short of the ones for wheat by 
around 20%. Maize at the same specific weight as wheat 
(Y = 0.80) produced a larger bottom pressure by 22% due to 
its smoother grain surface. For a silo cell that is supposed to 
contain maize, the stiffness of the cell walls and the bottom 
must be increased by 22%. 

The experiments described above - from which the cal- 
culations for arbitrary cell sizes are deduced - are performed 
in small apparatuses only, because the creation of larger sam- 
ple cells is connected with non-negligible costs. Regardless, 
they should be clarifying in some aspect, and it is hoped that 
they find confirmation in experiments on larger scale like the 
ones announced by a well-known mill-building firm. I shall 
follow up with a further communication on that subject at 
that point. 



2 Janssen and Predecessors 

Janssen cites only two earlier works, and these are given 
rather incompletely. However, they are believed to be cited 
correctly here. The first is Ref. [2], Die Silospeicher von 
Galatz und Braila - The silos of Galali and Brdila, where 
Galati and Braila are two Romanian cities at the Danube 
river. The second reference is [3], The Pressure of Stored 
Grain, and can be seen as an earlier account of the satura- 
tion effect observed in these silo experiments. This work is 
clearly acknowledged by Janssen, so it may seem unjustified 
to call this saturation effect alone the Janssen effect. 

There exists an earlier partial translation of Janssen's pa- 
per in [4], On the Pressure of Grain in Silos, which is listed 
there under the category of Foreign Abstracts and contains a 
short summary of the article by an editor. According to this 
editor, the saturation of the pressure "[...] has long been 
recognised". While Janssen's experiments are described in 
somewhat more detail, the Eqs. (2) and (3) are only cited 
without derivation and are indeed copied wrongly from the 
original. 

In fact, one can infer from earlier work that is not cited 
by Janssen, that the saturation of the pressure with depth was 
already more or less well known in Janssen's time. Janssen's 
work is anticipated in 1852 by a paper of G.H.L. Hagen [5], 
who became famous for his work in fluid dynamics, and 
an even earlier investigation by Huber-Burnand in 1829 [6]. 
Hagen imagined a cylinder of radius r and height h inside a 
container filled with sand of a specific weight y. The pressure 
exerted on the bottom area Kr 2 by the weight of the sand- 
cylinder is diminished by the friction experienced by that 
cylinder from the sand surrounding it. Assuming isotropic 
horizontal pressure throughout the system, Hagen derives 
the formula 

p(h) oc r 2 nyh - 2r%ylh 2 (4) 

with some friction coefficient I. He discards the decreasing 
branch of the solution and predicts the saturation of the pres- 
sure with depth after the height reaches the value h = r/(4l). 
Using the latter relation, Hagen subsequently tests his pre- 
dictions in an apparatus quite similar to Janssen's setup. He 
determines values for I between 0.154 and 0.22, but he does 
not check if the pressure below the maximum follows Eq. (4) 
which predicts a steeper increase than Janssen's Eq. (2). Ha- 
gen rather continues to analyze in detail the flow through 
an opening of his container and finds that the mass flow 
is proportional to some effective radius to the power 5/2, 
cf. Ref. [7, sec. 10.2]. For both flow and pressure measure- 
ments, Hagen mentions earlier work by Huber-Burnand with- 
out giving any reference. The appropriate reference is most 
likely Ref. [6] which discusses largely qualitative results of 
experiments with sand, peas, and - an egg. The egg is placed 
in a box, covered with several inches of sand, and loaded 
with a weight of 25 kg; the egg demonstrates - by not break- 
ing, of course - that only a small fraction of the pressure 
from the top actually reaches the bottom. 
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Fig. 15 Citations of the paper [1] since 1978 [9]. 



3 Impact of Janssen's paper 

The journal where Janssen published his results can still 
be found today [8] - it is still a publication of the Verein 
Deutscher Ingenieure - the Association of German Engi- 
neers. Figure 15 shows the citations of Janssen's paper [1] in 
the years since 1978. As of September 1st 2005, the article 
was cited 375 times with only 40 citations recorded for the 
time 1977 and earlier [9]. This demonstrates growing inter- 
est in granular systems in general and in one of its earliest 
papers in particular. 

It is interesting to note that Janssen's formula (2) is usu- 
ally derived for the more symmetric cylindrical geometry 
[7, sec. 5.2] and that more involved profiles require some 
averaging of the wall stress along the perimeter, which can 
be done explicitly, cf. [7, sec. 5.5], or implicitly as in [10, 
sec. 3.1.4]. In any case, the profile does only change the pref- 
actors in Eq. (2) but the overall behavior remains the same. 
Since he was experimenting with a square profile, Janssen 
was well aware of that kind of complication as seen in his 
more involved argument connected to Fig. 12. It was also 
found empirically that the saturation pressure was propor- 
tional to the diameter of the circle inscribed into the pro- 
file by Roberts [3]. Janssen reports this earlier finding, and 
Fig. 12 may suggest that he considered it essentially correct. 

Janssen's model rests on a number of approximations 
that were scrutinized later on: First, Eq. (1) amounts to full 
mobilization of friction at the walls. While that is not true 
in general, it might well hold to a good degree for Janssen's 
experimental protocol where the walls of the container were 
shifted upwards by the screws S in Figs. 1 and 2 after each 
step before filling in more grain. Second, the vertical and 
horizontal stresses are introduced as principal stresses which 
is clearly inconsistent with the existence of a finite friction 
at the walls. These principal stresses - horizontal and ver- 
tical pressures p s and p, respectively - are then related by 
a constitutive law, pK = p s f, which is not motivated any 
further. Third, a uniform distribution of stresses is assumed 
for any slice of the silo, cf. Fig. 10. Taking the constitutive 



law for granted, this translates into a finite shear stress in the 
center of the silo which must be ruled out by symmetry. All 
three of these objections can be overcome by more refined 
continuum theories [7], and the result in Eq. (2) remains 
valid with appropriate changes to the prefactors. There are, 
however, more serious problems when considering relatively 
large overloads on top of the grains in the silo [7, sec. 5.5]. 

If one probes in more detail the granular structure within 
the silo and especially when considering the dynamics of 
the granular material, Janssen's approach no longer yields 
satisfactory answers [10]. That said, it is quite remarkable 
that Janssen already postulates that it is arching within the 
granular assembly what prevented him from measuring the 
bottom pressure at various points. 
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